Divalent cations potentiate TRPV1 channel by lowering the heat 
activation threshold 
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Transient receptor potential vanilloid type 1 (TRPV1) channel responds to a wide spectrum of physical and 
chemical stimuli. In doing so, it serves as a polymodal cellular sensor for temperature change and pain. Many 
chemicals are known to strongly potentiate TRPV1 activation, though how this is achieved remains unclear. In 
this study we investigated the molecular mechanism underlying the gating effects of divalent cations Mg 2+ and 
Ba 2+ . Using a combination of fluorescence imaging and patch-clamp analysis, we found that these cations potenti- 
ate TRFVT gating by most likely promoting the heat activation process. Mg 2+ substantially lowers the activation 
threshold temperature; as a result, a significant fraction of channels are heat-activated at room temperature. 
Although Mg 2+ also potentiates capsaicin- and voltage-dependent activation, these processes were found either to 
be not required (in the case of capsaicin) or insufficient (in the case of voltage) to mediate the activating effect. 
In support of a selective effect on heat activation, Mg 2+ and Ba 2+ cause a Ca + -independent desensitization that 
specifically prevents heat-induced channel activation but does not prevent capsaicin-induced activation. These 
results can be satisfactorily explained within an allosteric gating framework in which divalent cations strongly pro- 
mote the heat-dependent conformational change or its coupling to channel activation, which is further coupled 
to the voltage- and capsaicin-dependent processes. 



INTRODUCTION 

Transient receptor potential vanilloid type 1 (TRPV1 ) is a 
heat- and ligand-sensitive ion channel that is activated 
when the temperature reaches ^40°C at resting mem- 
brane potential or when capsaicin is present at 100 nM 
or higher concentrations (Caterina et al., 1997). The 
dual activation by heat and capsaicin provides a molecu- 
lar explanation for the hot sensation elicited by chili 
peppers as well as the potential physiological role TRFVT 
plays in sensing temperature change and pain. In addi- 
tion, activity of TRPV1 is regulated by physical stimuli 
including voltage and mechanical force as well as a pleth- 
ora of chemicals such as extracellular H + , intracellular 
Ca 2+ , and PIP 2 (Tominaga et al., 1998; Jordt et al., 2000; 
Chuang et al., 2001; Voets et al., 2004; Stein et al., 2006; 
Lukacs et al., 2007; Dhaka et al., 2009; Ufret-Vincenty 
et al., 2011; Cao et al., 2013). Sensitivity to a wide spec- 
trum of stimuli allows TRPV1 to serve as a polymodal 
cellular sensor (Clapham, 2003; Zheng, 2013). Under- 
standing how TRPV1 senses these stimuli is of great 
practical importance, as the channel is considered to 
be an attractive drug target for pain medication (Wu 
etal., 2010). 
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Like voltage-gated potassium channels, TRPV1 is 
a tetrameric protein complex with a centrally located 
ion permeation pore surrounded by channel subunits 
that contain six transmembrane segments and intracellu- 
larly located amino and carboxyl termini (Jahnel et al., 
2001; Kedei et al., 2001; Kuzhikandathil et al., 2001; 
Moiseenkova-Bell et al., 2008). Chemical activators of 
TRFVT are known to interact with many different channel 
structures. Capsaicin, for example, interacts with resi- 
dues in the S2-to-S3 region, including likely tt-tt inter- 
actions between the vanillyl moiety of capsaicin and the 
aromatic ring of amino acids Y511 and/or F512 (of rat 
TRPV1; Jordt and Julius, 2002). Ca 2+ -calmodulin binds 
to two potential intracellular binding sites, located in 
the C-terminal region and the N-terminal Ankyrin-like 
repeat domain (Lishko et al., 2007; Lau et al., 2012). 
Extracellular H + is found to bind to sites clustered at the 
outer pore region, mainly two glutamates (E600 and 
E648 in rat TRPV1; Jordt et al., 2000) . The wide spread 
of chemical interaction sites suggests that an allosteric 
mechanism may underlie the integration of multiple 
TRPV1 stimuli in promoting channel activation. Gating 
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models incorporating allosteric mechanisms indeed 
have been proven to be successful in describing many 
aspects of TRPV1 gating (Latorre et al., 2007; Matta and 
Ahern, 2007; Jara-Oseguera and Mas, 2013). 

In contrast to the detailed understanding of chemical 
activation of TRPV1, the molecular mechanisms gov- 
erning the channel's response to most physical stimuli 
remain unclear. In particular, the heat activation pro- 
cess that underlines the role of TRPV1 as a heat and 
pain sensor is highly controversial. An important factor 
impeding research in this area is the lack of effective 
methods to pinpoint the sites affected by heat. It is well 
known that structural perturbations by mutation (dele- 
tion, chimera, etc.) of selected channel regions do not 
have only localized effects. Furthermore, unlike chemi- 
cal stimuli that target specific protein sites, stimulus 
from heat is difficult to restrict to a specific protein site 
or region. When temperature rises, heat affects the 
whole protein as well as its surrounding lipids and the 
aqueous environment. As channel activity can be poten- 
tially affected by each of these components, it is difficult 
to distinguish what changes drive channel activation. 

Given the available knowledge on the action sites and 
molecular mechanisms for TRPV1 chemical stimuli, it 
would be informative to know how chemical stimuli in- 
teract with the physical stimuli of TRPV1 , especially if 
such interaction is specific. Results from our recent stud- 
ies (Yang et al., 2010; Cui et al., 2012) and other studies 
(Matta and Ahern, 2007; Grandl et al., 2010) revealed 
that capsaicin and heat apparently activate TRPV1 
through separate pathways. In this paper, we investi- 
gated the gating mechanism of another type of chemi- 
cal stimuli for TRPV1 , the divalent cations. Divalent and 
multivalent cations such as Mg 2+ , Ni 2+ , Gd 3+ , and poly- 
amines are potent TRPV1 activators (Ahern et al., 2005, 
2006; Tousova et al., 2005; Riera et al., 2007; Ohta et al., 
2008; Luebbert et al., 2010). Divalent cations have sim- 
ple chemistry, high water solubility, and low membrane 
permeability (in contrast to, for example, capsaicin), 
and are chemically inert. These preferable features make 
divalent cations an attractive tool for the investigation 
of the TRPV1 gating mechanism. For these reasons, in 
this paper and the accompanying paper (see Yang et al. 
in this issue) we examined the mechanism and structural 
determinants of TRPV1 potentiation by two divalent 
cations, Mg 2+ and Ba 2+ . 

The present study focused on the activation mecha- 
nism. We assessed the permeation effect of Mg 2+ , the 
concentration dependence in promoting channel acti- 
vation, and the channel open probability that could be 
reached by Mg 2+ stimulation. We then quantified and 
compared the energetic effects of Mg 2+ on various TRPV1 
activation pathways. Our observations showed that these 
cations strongly potentiate channel activation, which 
can be satisfactorily explained within an allosteric frame- 
work assuming that divalent cations directly potentiate 



the heat-activation pathway, but cannot be easily ex- 
plained if instead their gating effects are assumed to 
originate from potentiating the voltage- or capsaicin- 
activation pathways. In support of a specific effect on 
heat activation, we observed that prolonged application 
of Mg 2+ and Ba 2+ selectively desensitizes heat activation 
while sparing capsaicin activation. These results identify 
divalent cations as a distinct group of TRPV1 agonists. 
In addition, they further support the notion that heat 
and capsaicin activate TRPV1 through different pathways 
(Jordt et al., 2000; Matta and Ahern, 2007; Yang et al., 
2010; Cui et al., 2012). Our findings demonstrate that 
divalent cations can be a useful tool in the search for the 
heat activation mechanism. In Yang et al. (2014) , we de- 
scribe experiments aimed at identifying channel struc- 
tures that mediate the gating effects of divalent cations. 

MATERIALS AND METHODS 

cDNA constructs and cell transfection 

The mouse TRPV1 (a gift from M.X. Zhu, University of Texas 
Health Science Center at Houston, Houston, TX) was used in 
this study. This cDNA has 95% and 87% overall sequence identity 
to the rat and human TRPVl cDNAs, respectively. Functional 
properties of mTRPVl are also very similar to rTRPVl and hTRFVT, 
on which several previous studies on the gating effects of cationic 
ions were performed. To facilitate identification of channel- 
expressing cells, the cDNA encoding enhanced YFP (eYFP) was 
fused in frame to the C-terminal end of the mTRPVl cDNA, as 
described previously (Cheng et al., 2007). 

HEK293 cells were cultured in a DMEM medium supplemented 
with 10% FBS, 100 U/ml penicillin, and 100 mg/ml streptomycin 
at 37°C with 5% C0 2 . Cells were passaged 6-24 h before trans- 
fection by plating onto glass coverslips coated with 0.1 mg/ml 
poly-L-lysine to improve cell adhesion and subsequent patch 
recordings. Transient transfection was conducted by adding 
4 pg plasmid DNA and 4 pi Lipofectamine 2000 (Invitrogen) 
in a 6-well plate. Electrophysiological experiments and fluores- 
cence imaging recordings were performed between 24 h and 
48 h after transfection. 

Electrophysiological recordings 

Macroscopic and single-channel currents were recorded from 
TRPVl-expressing cells using a HEKA EPC10 amplifier controlled 
with PatchMaster software (HEKA). Patch pipettes were pulled 
from borosilicate glass and fire-polished to a resistance of ~2 Mfi. 
Membrane potential was held at 0 mV and, unless otherwise 
stated, currents were elicited by a protocol consisting of a 300-ms 
step to +80 mV followed by a 300-ms step to —80 mV at 1-s inter- 
vals. Data were filtered at 2.25 kHz and sampled at 12.5 kHz. As 
TRPVl macroscopic current was strongly outward-rectifying mainly 
due to rapid deactivation at — 80 mV, current amplitude was ana- 
lyzed at +80 mV. Both whole-cell and inside-out configurations 
were used in this study, while Mg 2+ or Ba 2+ were always added to 
the extracellular solution. During whole-cell recordings the ca- 
pacitance current was minimized by amplifier circuitry. The series 
resistance (typically 5-20 Mf2) was compensated by 80%; after 
compensation the voltage errors should be <5 mV. The standard 
bath solution and pipette solution contained 130 mM NaCl and 
3 mM Hepes, pH 7.2. Mg 2+ or Ba 2+ solution contained 130 mM 
MgCl 2 or BaCl 2 , and 3 mM Hepes, pH 7.2. Solutions containing 
10-100 mM Mg 24 were prepared from the 130 mM Mg 2+ -containing 
solution by replacing the appropriate concentration of MgCl 2 
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with NaCl. Solution switching was achieved with a rapid solution 
changer RSC-200 (Bio-Logic USA). 

Temperature control and monitoring 

The bath solution was heated using an SHM-828 eight-line heater 
driven by a CL-100 temperature controller (Harvard Apparatus). 
To obtain a complete temperature-current relationship, in some 
experiments the bath solution was first cooled by embedding the 
perfusion solution reservoir in ice water. We placed a TA-29 min- 
iature bead thermistor (Harvard Apparatus) ~1 mm from the 
pipette tip to monitor local temperature change. Temperature 
readout from the thermistor was fed into an analogue input port 
of the patch-clamp amplifier and recorded simultaneously with 
channel current. The speed of temperature change was set at a 
moderate rate of ^0.3°C/ s to ensure that heat activation reached 
equilibrium during the course of temperature change and the 
current was recorded at steady-state. When the experimental tem- 
perature was not controlled, recordings were conducted at room 
temperature at 24°C. 

Using fluorescence imaging to record TRPV1 activity 

To monitor TRPV1 activity in intact living cells, we chose to use 
the high-affinity calcium indicator Fluo-4 (i^ for Ca 2+ is 345 nM), 
which yields a large fluorescence intensity increase upon Ca 2+ 
binding and its excitation spectrum, peaking at 494 nm, and is 
compatible with an argon laser (Gee et al., 2000). As a divalent 
cation chelator, Fluo-4 may potentially also bind Mg 2+ , especially 
as in our experiments concentrations up to 130 mM (comparing 
to 1.8 mM Ca 2+ ) were applied to cells expressing TRPV1, a nonse- 
lective cation channel. Although there is no report available on 
the binding affinity of Mg 2+ to Fluo-4, Mg 2+ binds to the closely 
related Ca 2+ indicator dye Fluo-3 with an extremely low affinity 
(Ka of 9 mM; Minta et al., 1989) . This low binding affinity, in com- 
bination with the observation that binding of Mg 2+ to Fluo-3 in- 
creases fluorescence weakly (only 1.4-fold, comparing to a 40-fold 
increase by Ca 2+ binding), suggests that in our experiments the 
fluorescence signal came predominantly from binding of Ca 2+ 
that entered the cell through open TRFVT channels, whereas 
Mg 2+ might also contribute to the total fluorescence change. This 
prediction was confirmed experimentally. As shown in Fig. SI, in 
the absence of Ca 2+ , Mg 2+ was indeed able to yield small fluores- 
cence signals when entering the cell, but only <10% in amplitude 
comparing to the signal generated by Ca 2+ . Control experiments 
using cells pretreated with thapsigargin (1 uM), an inhibitor for 
sarcoplasmic reticulum Ca 2+ ATPase (SERCA pump), did not ex- 
hibit a noticeable difference in fluorescence signal intensity or 
kinetics compared with untreated cells. This observation further 
confirmed that changes in fluorescence intensity were caused by 
Ca 2+ entry upon heat, capsaicin, or Mg 2+ stimuli, instead of the 
release of Ca 2+ from the intracellular store. 

For fluorescence imaging experiments, TRPVl-expressing 
HEK293 cells seeded on 25-mm coverslips were subjected to 
imaging 24-48 h after transfection. After washing once with an 
extracellular solution (ECS) containing 140 mM NaCl, 5 mM KC1, 

1 mM MgCl 2 , 1.8 mM CaCl 2 , 10 mM glucose, and 15 mM HEPES, 
pH 7.4, cells were incubated in 2 ml of ECS supplemented with 

2 uM fluo-4/ AM and 0.1% Pluronic F-127 (both from Invitrogen) 
at room temperature for 60 min. 2 mM Probenecid (Invitrogen) 
was included in all solutions during imaging to prevent fluo-4 
leakage from loaded cells. After incubation, cells were washed 
three times with ECS and incubated in the same solution for an- 
other 20 min to complete the intracellular hydrolysis process of 
the AM ester, which converts the nonfluorescent Fluo-4/AM into 
the fluorescent Fluo-4. 

Coverslips with dye-loaded cells were placed in the quick-re- 
lease magnetic chamber (Warner Instruments) and mounted 
on the stage of a microscope system (Eclipse TE2000-U; Nikon) 



equipped with a CCD camera (Cascade 128B; Roper Scientific). 
Fluo-4 was excited by an Argon laser with a filter set of Z488/10 
(excitation), z488rdc (dichroic) and recorded through an emis- 
sion filter HQ5001p (all from Chroma Technology Corp.). The 
same excitation intensity was used for all the experiments. The 
duration of light exposure was controlled by a computer-driven 
mechanical shutter (Uniblitz). Cell images were acquired sequen- 
tially with an exposure period of 200 ms at an interval of 1 s. Both 
the shutter and the camera were controlled with MetaMorph soft- 
ware (Universal Imaging) . 

ECS with varying concentration of Mg 2+ was delivered right 
next to the imaging area by an RSC-200 rapid solution-switching 
system through separate lines. Temperature control during imag- 
ing was achieved by perfusing preheated solutions in the same 
way as for electrophysiology studies. A TA-29 miniature bead therm- 
istor (Harvard Apparatus) was placed right next to the imaging 
area to ensure accurate monitoring of local temperature. The 
thermistor's temperature readout was fed into an analogue input 
of the patch amplifier; MetaMorph software was synchronized 
with PatchMaster software to record temperature simultaneously 
with imaging. All experiments except heat activation studies were 
performed at room temperature. 

For imaging data processing, the cell boundary was defined 
manually in MetaMorph. The boundary defined individual re- 
gions of interest (ROIs). The mean fluorescence intensity of each 
ROI was measured and exported to Excel (Microsoft) . When an 
equilibrium state was reached after application of a stimulus, the 
mean fluorescence intensity, F, of all pixels within each ROI was 
calculated. Changes in fluorescence intensity, AF, were calculated 
as the difference between the equilibrium level after stimulus ap- 
plication and the baseline. By measuring fluorescence changes, 
background fluorescence from eYFP and scatter light were ex- 
cluded. The background fluorescence was not affected by Mg 2+ 
under our experimental conditions. 

Data analysis 

We used the current-temperature relationship to determine acti- 
vation threshold temperature and to characterize thermodynamic 
properties of heat activation (Cui et al., 2012). The raising phase 
of the current-temperature curve recorded from cells expressing 
wild-type or mutant TRPV1 channels showed two temperature- 
dependent phases: a less temperature-dependent phase at lower 
temperatures followed by a higher temperature-dependent phase 
at higher temperatures. We fitted each phase to a linear function. 
The temperature at the intersection of the two lines was defined 
as the activation threshold temperature. The same approach was 
used to estimate the activation threshold from fluorescence imag- 
ing data. 

To calculate the enthalpic change, AH, and the entropic 
change, A S, associated with heat activation, we chose the high tem- 
perature-dependent phase to construct a Van't Hoff plot and fit- 
ted it to the following equation: 

AH AS 

mK = + — , 

q RT R 

where K eq is the equilibrium constant calculated from the channel 
open probability, R is the gas constant, and T is the temperature 
in Kelvin. Because the single-channel conductance is tempera- 
ture dependent, the macroscopic current amplitudes at different 
temperatures were corrected by converting them to a level as if 
the temperature were 20°C (Yang et al., 2010). The open proba- 
bility was determined as the ratio between the macroscopic cur- 
rent (after correcting for temperature-dependent single-channel 
conductance) and the maximum current estimated using noise 
analysis (Sigworth, 1980). The mean current amplitude, /, and 
the corresponding variance, cr 2 , were estimated from the plateau 
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phase of each sweep. Noise analysis was then applied to estimate 
the number of channels and the maximum current level, using 
the following equation o~ 2 = il • m which i is the single- 

channel amplitude and N is the number of channels. The half- 
activation temperature, T^ip at which the heat activation is 50% 
complete, was calculated as T llalf = ■ 

G-Vcurves were fitted to a single-Boltzmann function: 



1 + e 



1 



(1) 



where G/G mrlx is the normalized conductance estimated from the 
steady-state current at each testing voltage, q is the equivalent gating 
charge, V^^is the half-activation voltage, and iris Faraday's constant. 

The capsaicin dose-response relationship was fitted with the 
Hill equation: 



[xl 



EC 5 V[x] n 



where I x is the steady-state TRPV1 current in the presence of cap- 
saicin at concentration [x] , I m i„ and I max are the current amplitude 
in the absence and the presence of a saturating concentration of 
capsaicin, respectively, EC; 0 is the capsaicin concentration at 
which activation is half-maximal, and n is the slope factor. 

Single-channel conductance was estimated from all-point his- 
tograms constructed from current events recorded at +80 mV. 
A double-Gaussian function was fitted to the histograms. From 
the positions of the Gaussian peaks the closed and open current 
levels were identified, and single-channel current amplitude and 
the corresponding conductance were calculated. 

All values are given as mean ± SEM for the number of measure- 
ments indicated ( n) . Statistical significance was determined using 
the Student's t test, and indicated as follows: *, P < 0.05; **, P < 
0.01; ***,P< 0.001. 

Multi-allosteric gating model 

The approach we used for the study of the TRPV1 polymodal gat- 
ing was inspired by the work on Ca 2+ -sensitive voltage-gated BK 
channels by Aldrich, Horrigan, and Cui (Horrigan and Aldrich, 
1999, 2002; Horrigan et al., 1999). Using a structure-indepen- 
dent, dual-allosteric gating model, these authors demonstrated 
nicely that near-independent gating transitions caused by Ca 2+ 
binding and voltage sensor movement exert substantial influence 
on each other through their convergence onto the channel open- 
ing transition. Similarly, our approach to a general gating model 
for TRPV1 assumed that the closed-to-open transition, C<->0, and 
the transition induced by a given activation stimulus, R<->A, are 
both reversible transitions that are coupled allosterically (Scheme 
I in Fig. 1). In this case, the channel open probability, P„, is deter- 
mined by the equilibrium constants of the two transitions, L and 
D, as well as the allosteric coupling factor between the two transi- 
tions,/, according to the following equation: 



m+jp) 

L(\ + JD) + (l + D)' 



(2) 



where L = io exp(wi/i?7) and D = D 0 exp(w D /RT). w L and w D sxe 
the free energy associated with the C<->0 and R<-»A transition, 
respectively. Through allosteric coupling factor J, the equilib- 
rium constant D exerts its effect on the channel open probability, 
which is the experimentally observable parameter. For a strong 
coupling, J S> 1, and P 0 is strongly dependent on D. In the ex- 
treme case when J = 1, there is no coupling between the two 
transitions, and the open probability becomes independent of 
D-.P 0 =L/(1+L). 



Applying this general, structure-independent principle to the 
polymodal activation of TRPV1 by voltage, capsaicin, and heat 
yielded Scheme II (Fig. 1 B), in which the closed-to-open transi- 
tion, GhO, is allosterically coupled to the voltage-dependent 
transition between the resting state and the activated state, R<->A, 
the capsaicin-induced transition between the unliganded state 
and the liganded state, Uf>L, and the heat-induced transition be- 
tween the native state and the permissive state, N<->P, with a cou- 
pling factor of Jv, Jc, and Jh, respectively. When two instead of 
three activation stimuli are present simultaneously, Scheme II 
(Fig. 1 B) reduces to a dual-allosteric system equivalent to the one 
for the BK channel (Horrigan and Aldrich, 2002) . Under condi- 
tions that prohibit all but one activation pathway, Scheme II re- 
duces to Scheme I (Fig. 1), with the channel open probability 
determined by two equilibrium constants and one coupling factor 
as shown by Eq. 2. In support of the multi-allosteric gating model 
of Scheme II (Fig. 1 B), we and others have found that heat and 
capsaicin activate TRPV1 through separate pathways that are en- 
ergetically coupled (Grandl et al., 2010; Yang et al., 2010; Cui 
et al., 2012). Similarly, voltage and capsaicin stimuli exhibit an 
additive effect on activation (Matta and Ahern, 2007; Yang et al., 
2010). In principle, transitions induced by the three stimuli can 
also be allosterically coupled (as shown by the dashed lines in Fig. 1 
B, Scheme II) with coupling factors Jcujcv, and Jm- Observations 



A Scheme I. allosteric coupling 



R 

laJ 



B Scheme II. TRPV1 gating pyramid 
N 




Figure 1. Allosteric gating schemes. (A) A general allosteric 
coupling system containing a reversible closed (C) to open (O) 
transition with an equilibrium constant L, a reversible transition 
between a relaxed state (R) and an activated state (A) with an 
equilibrium constant D, and a coupling factor J. The measurable 
parameter is the probability of being in the O state, which is deter- 
mined by both the C<->0 equilibrium and the R<->A equilibrium. 
No additional constraint is assumed in this structure-independent 
model. (B) A multi-allosteric model for the gating of TRPV1 by 
voltage, capsaicin, and heat, in which the C<->0 equilibrium at 
the apex of the pyramid is coupled allosterically to each of the 
stimulus-induced transitions according to the general principle 
shown in Scheme I. Definition of states and coupling factors are 
given in Materials and methods. 
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from our previous experiments suggested that Jch is apparently 
quite small (Cui etal., 2012). Nonetheless, in this study we did not 
attempt to estimate the strength of coupling between different 
stimuli. Likewise, we were not concerned with obtaining an accu- 
rate estimate of the coupling factors J v , J c , and J H . Our focus was 
instead on estimating the energetic effect of divalent cation Mg 2+ 
on each of the three branches of Scheme II (Fig. 1 B) , from which 
we can gain insights into the mechanistic nature of Mg 2+ potentia- 
tion of TRPV1. The energetic effect of 100 mM Mg 2+ on the volt- 
age-dependent activation was estimated from the shift of the G-V 
curve upon application of the divalent cation: AE V = qFAV, in 
which AEy is the free energy difference in voltage-dependent acti- 
vation in the absence and presence of Mg 2+ , qis the apparent gating 
charge, F is the Faraday's constant, and A V is the shift in half- 
activation voltage Vkaif (see Eq. 1). The energetic effect of 10 mM 
Mg 2+ on the capsaicin-induced activation was estimated from the 
shift of the capsaicin dose-response curve: AE C = 7iRTln(EC 5Q / 
EC m) , in which A£ c is the free energy difference in ligand activa- 
tion caused by Mg 2+ , EQ,o and EC 50 are the half-activation concen- 
trations of capsaicin in the absence and presence of Mg 2+ , 
respectively, and n is the Hill slope factor. The energetic effect of 
10 mM Mg 2+ on heat-dependent activation was estimated from the 
change in the free energy associated with the heat activation pro- 
cess: AE H = A — AGi, in which AE H is the free energy difference 
in heat activation in the absence and presence of Mg 2+ , AGi, and 
A&. AGi and AGj were calculated from the enthalpy and entropy 
using the Gibbs equation A G = AH — TAS as described earlier. 

Online supplemental material 

Fig. SI shows that fluorescence intensity increase upon channel 
opening was mainly induced by Ca 2+ , whereas Mg 2+ could also 
induce a small fluorescence signal. Fig. S2 shows that the rep- 
resentative time course of the temperature-dependent fluores- 
cence intensity changes upon heating in a TRPVTexpressing 
cell. Online supplemental material is available at http://www.jgp 
.org/cgi/content/full/jgp.201311025/DCl. 



RESULTS 

Mg 2+ induces influx of divalent cations in 
TRPV1 -expressing HEK293 cells 

We first examined the effect of Mg 2+ on TRFV1 activa- 
tion by recording intracellular fluorescence increase 
due to Ca 2+ (and Mg 2+ ) influx through open channels. 
The mouse TRPV1 channel was expressed in HEK293 
cells that were then loaded with the Ca 2+ indicator dye 
Fluo-4. Introducing Mg 2+ into the bath solution at room 
temperature induced a fluorescence signal in a dose- 
dependent manner (Fig. 2 A) . Although 30 mM Mg 2+ 
could induce a clearly observable increase in fluores- 
cence intensity, further increasing the Mg 2+ concentra- 
tion to 100 mM significantly boosted the fluorescence 
signal, indicating that the apparent EC50 value must be 
>30 mM (Fig. 2, B and C). It is worth noting that this is 
at least three-to-four orders of magnitude higher than 
the EC50 value for capsaicin on mouse TRPV1, which we 
previously estimated to be 0.95 pM (Cui et al., 2012), as 
well as the EC50 value for proton, at ^10 pM, pH 5 (un- 
published data) . The Mg 2+ sensitivity of mouse TRPV1 
appeared to be slightly lower than that of human 
TRPV1 recorded with a different Ca 2+ dye, Fura-2 (Riera 



et al., 2007). We noticed that, with 100 mM Mg 2+ , the 
fluorescence intensity declined over time with a time 
constant of many seconds (see Fig. 2 B), which is cor- 
related with a slow desensitization process (as described 
later) . In summary, our results from intact live cells 
suggest that at millimolar concentrations, Mg 2+ causes 
TRPV1 activation at room temperature, followed by 
slow desensitization. 

Mg 2+ strongly potentiates TRPV1 current 

To directly observe Mg 2+ effects on TRFVT activation, 
we recorded TRFVT currents using whole-cell patch 
clamp in the presence of various concentrations of 
Mg 2+ in the bath. Channel-expressing cells were iden- 
tified by the fluorescence signal from an eYFP tag at 
the channel's C-terminal end, which were confirmed 
by the elicitation of TRPVl-mediated currents by 3 pM 
capsaicin. In close agreement with previous reports on 
rat and human TRPV1 channels (Ahern et al., 2005; 
Wang et al., 2010), extracellular Mg 2+ at 10-100 mM 
concentrations strongly potentiated TRPV1, eliciting 
noticeable currents at room temperature (Fig. 3 A). 
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Figure 2. Mg 2+ induces Ca 2+ influx in TRPVTexpressing cells. 
(A) Live-cell fluorescence imaging of changes in intracellular 
Ca 2+ concentration. (B) Representative time courses of the fluo- 
rescence signal in a channel-expressing cell and an untransfected 
cell. Changes in fluorescence intensity are shown in arbitrary 
units after subtracting the basal level intensity. (C) Change in 
fluorescence intensity increases with Mg 2+ concentration. ***, P < 
0.001; n = 70-1 13. Error bars indicate +SEM. 
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Consistent with fluorescence imaging recordings, cur- 
rent induction by Mg 2+ exhibited clear concentration 
dependence and did not show any sign of saturation at 
100 mM. At this concentration, the Mg 2+ -eli cited cur- 
rent was ~15% in amplitude compared with that elicited 
by 3 pM capsaicin, a saturating concentration (Fig. 3 B, 
black bars) . 

Mg 2+ reduces TRPV1 single-channel conductance 
We conducted single-channel recording in the inside- 
out configuration, with various concentrations of Mg 2+ 
in the pipette solution. We observed that extracellular 
Mg 2+ clearly made the channel opening events longer 
and occur more frequendy (Fig. 4 A) , which is consis- 
tent with a potentiating effect on gating. In addition, 
Mg 2+ dose-dependently inhibited the single-channel 
current amplitude (Fig. 4, A and B; Samways and Egan, 
201 1) . In the Mg 2+ -free solution, the single-channel con- 
ductance of TRPV1 outward current was estimated to be 
133.8 ± 3.4 pS (n = 4), similar to that reported previ- 
ously (Caterina et al., 1997; Premkumar and Ahern, 
2000) . At the 100-mM concentration, Mg 2+ reduced sin- 
gle-channel conductance by 50%, to 66.5 ± 0.1 pS (n = 4; 
Fig. 4, B and C). Thus Mg 2+ dose-dependently potenti- 
ates TRPV1 through a gating effect but at the same time 
inhibits conductance. The dual effect of Mg 2+ on gating 
and permeation could be seen in a previous study in 
which Mg 2+ potentiated human TRPV1 current elicited 
by capsaicin (which partially activated the channel) but 
inhibited rat TRFVT current elicited by capsaicin (which 
fully activated the channel; Wang et al., 2010) . 

Estimation of the channel open probability 
induced by Mg 2+ 

As the permeation effect of Mg 2+ would partially mask 
the gating effect, in order to estimate the open proba- 
bility of TRPV1 in the presence of Mg 2+ we used the 
Mg 2+ concentration dependence in single-channel con- 
ductance to adjust the whole-cell current amplitudes. 
Removing the permeation effect revealed a larger gat- 
ing effect by Mg 2+ . We estimated that, at saturating cap- 
saicin concentrations and +80 mV, the mouse TRFVT 
activation reaches an open probability of ^80%. Based 
on the amplitude ratio between the adjusted Mg 2+ -elicited 




current and the capsaicin-elicited current, 100 mM Mg + 
potentiated the channel to an open probability of ^25% 
(n = 5; Fig. 3 B, red bars). Our data thus demonstrated 
that Mg 2+ is a strong activator for TRFV1 . 

Having established that Mg 2+ directly and potently 
activates TRPV1, we focused our attention on the mo- 
lecular mechanism underlying the Mg 2+ -mediated acti- 
vation process. Our approach was based on a general, 
structure-independent, multi-allosteric framework in 
which each of the major TRFV1 stimuli (capsaicin, volt- 
age, and heat) promotes a transition in the channel that 
is allosterically coupled to the closed-to-open transition 
(Fig. 1 B) . We wanted to test which activation pathway 
is predominantly affected by Mg 2+ . For this purpose we 
designed experiments to examine the energetic effect 
of Mg 2+ on each of the activation pathways. 

Effects of Mg 2+ on voltage-dependent activation 
One potential mechanism for Mg 2+ potentiation of 
TRFV1 gating is through shifting the voltage-dependent 
activation process. To investigate this possibility, we com- 
pared the voltage dependence of TRFVT activation with 
and without extracellular Mg 2+ (Fig. 5 A). At 100 mM, 
Mg 2+ shifted the G-V curve by —46 mV without a sig- 
nificant change in the steepness (Fig. 5 B). The shift in 
Vhaif would make depolarization more effective in acti- 
vating TRPV1 . However, because TRPV1 is very weakly 
voltage-dependent, having an apparent gating charge 
of ^0.5 e 0 (Fig. 5 B; Voets et al., 2004; Matta and Ahern, 
2007), it is unlikely that the Mg 2+ potentiating effect is 
wholly caused by the shift in voltage dependence. In- 
deed, the shift in the G-V curve alone will increase the 
open probability at +80 mV by a mere 61% (Fig. 5 B), 
whereas based on the adjusted current amplitudes 
shown in Fig. 3 B, at +80 mV, 100 mM Mg 2+ increased 
the open probability by 26-fold. From the shift of V^aif 
and the overall voltage dependence of channel activa- 
tion, it was estimated that the energetic effect of 100 mM 
Mg 2+ on voltage-dependent activation is ^0.5 kcal/ mol. 
(The effect of 10 mM Mg 2+ would be too small to esti- 
mate accurately.) 
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Figure 3. Mg 2+ activates TRPV1 current in a dose- 
dependent manner. (A) Representative whole-cell 
currents at +80 mV evoked by Mg 2+ and capsaicin. 
(B) Mean currents (black bars) evoked by Mg 2+ 
normalized to the response evoked by 3 uM cap- 
saicin (left axis) . The red bars represent the same 
data after correction for Mg 2+ inhibition of the sin- 
gle-channel conductance (see Fig. 4), so that the 
height of the bars directly reflects relative open 
probability (right axis), n = 5 each. Error bars in- 
dicate ±SEM. 
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Figure 4. Mg 2+ reduces the single-channel conductance of TRPV1. (A) Representative single-channel current traces recorded at 
+80 mV with the indicated concentrations of extracellular Mg 2+ . (B) All-point histograms of single-channel events at the indicated Mg 2+ 
concentrations. The superimposed curve represents a fit of a double-Gaussian function. (C) Box-and-whisker plot of the single-channel 
conductance versus the corresponding Mg 2+ concentration. The whisker top, box top, line inside the box, box bottom, and whisker 
bottom represent the maximum, 75th percentile, median, 25th percentile, and minimum value of each pool of conductance measure- 
ments, respectively, n = 3-4. 



Effects of Mg 2+ on capsaicin activation 
Potentiation of TRPV1 by Mg 2+ does not require the 
presence of capsaicin. Nonetheless, Mg 2+ and capsaicin 
do exhibit synergy in activating the channel so that in 
the presence of Mg 2+ , capsaicin activation is more promi- 
nent (Ahern etal., 2005; Wang et al., 2010). We observed 
that Mg 2+ shifted the capsaicin dose-response curve ap- 
preciably to the left (Fig. 6), which would make cap- 
saicin a more potent agonist in the presence of Mg 2 *. 
To further characterize interaction between Mg 2+ and 
the capsaicin activation pathway, we estimated the en- 
ergetic effect of Mg 2+ on capsaicin activation from the 



Mg -induced change in ECso- As shown in Fig. 6, 10 mM 
Mg 2+ shifted the EC 50 value from 977 + 15 nM (n= 3) 
to 282 ± 20 nM (n = 4) with a minor effect on the Hill 
slope factor (1.88 ± 0.05 and 1.61 ± 0.10, respectively; P > 
0.05) . The shift in EC50 value reflects a 1.2 kcal/mol en- 
ergetic contribution to activation at room temperature. 

Mg 2+ substantially lowers the activation 
threshold temperature 

We next examined the effect of Mg 2+ on heat activation. 
Mg 2+ of various concentrations was added to the pi- 
pette solution, and heat-induced TRPV1 activation was 



Figure 5. Effects of Mg 2 * on the volt- 
age-dependent activation of TRFVT. 
(A) Representative whole-cell patch- 
clamp current traces in response to a 
family of voltage steps in the absence 
(left) or presence (right) of 100 mM 
Mg 2+ . Voltage steps were applied from 
a holding potential of 0 mV to vari- 
ous membrane potentials from — 1 60 
to +160 mV, in 20-mV steps. (B) G-V 
curves in the absence (open symbols) 
or presence (closed symbols) of 100 
mM Mg 2+ normalized to the respective 
G mo!t value. Superimposed are fits of a Boltzmann function with the following V M /and q values: 0 Mg 2+ , 130.5 ±1.6 mV, 0.53 ± 0.03 e 0 
(n = 5); 100 mM Mg 2+ , 84.3 ± 1.9 mV, 0.56 ± 0.06 e 0 (n = 5). Asterisks indicate a significant difference in the normalized conductance 
values (P < 0.05). Error bars indicate ±SEM. 
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Figure 6. Effects of Mg" + on the capsaicin-dependent activation of TRPV1. (A) Representative whole-cell current traces in response to 
varying capsaicin concentration in the absence (left) or presence (right) of 10 mM Mg 2+ . (B) Dose-dependent activation of TRPV1 by 
capsaicin in the absence (closed circles) or presence (open circles) of 10 mM Mg 2+ . Curves represent fits of a Hill equation. The EC;o 
and slope factor values are: 977 ± 15 nM and 1.88 ± 0.05 (n = 3) for 0 mM Mg 2+ , and 282 ± 20 nM and 1.61 ± 0.10 (n = 4) for 10 mM Mg 2+ . 
Error bars indicate ±SEM. 



recorded from inside-out patches perfused by a heated 
solution. Heat-induced channel activation was moni- 
tored by the time course of current increase. Unlike a 
previous observation with rat TRPV1 in oocytes (Ahern 
et al., 2005), Mg 2+ was found to substantially shift the 
temperature dependence of TRPV1 (Fig. 7, A and B) . 
In the absence of divalent cation, TRPV1 activated at 
37.7 ± 0.3°C (n = 9; Fig. 7, B and C) . With 10 mM Mg 2+ , 
the threshold temperature was significandy lowered to 
30.2 ± 1.1 °C (n = 8, P < 0.001). Increasing the Mg 2+ con- 
centration to 30 mM further shifted the threshold tem- 
perature to 20.2 ± 1.8°C (n = 8, P < 0.001) so that even 
at room temperature TRPV1 was open at an apprecia- 
ble level (Fig. 7 C). The large shift in threshold tem- 
perature made it impractical to further increase Mg 2+ 



concentration for these experiments. Nonetheless, it is 
clear that Mg 2+ has a strong effect on heat activation. 

Fluorescence imaging experiments confirmed that 
Mg 2+ substantially affected the heat response of TRPV1. 
While raising temperature in the presence of 1 mM 
Mg 2+ (the control solution) caused a gradual increase in 
the fluorescence signal, the same manipulation in the 
presence of 10 mM Mg 2+ elevated the fluorescence sig- 
nal to a higher level more rapidly and at lower tempera- 
tures (Fig. 8 A; see also Fig. S2). With 10 mM Mg 2+ , the 
peak fluorescence intensity was more than doubled 
(Fig. 8 B) , whereas the mean activation threshold tem- 
perature for fluorescence increase was shifted by M3°C 
(Fig. 8 C). These results from live cells again suggest 
that further increases in the Mg 2+ concentration would 
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Figure 7. Mg 2+ strongly potentiates the tempera- 
ture response of TRPV1. (A) Normalized current 
(top) induced by a temperature ramp (bottom) 
recorded from inside-out patches at +80 mV with 
0 mM, 10 mM, and 30 mM Mg 2+ in the pipette solu- 
tion. Recordings are from three separate patches. 
(B) The current raising phase is shown as a function 
of temperature. The activation threshold tempera- 
tures, at the intersection of a pair of dotted lines 
(representing the leak current and the heat-activated 
TRPV1 current), are indicated by red arrows. (C) Mg 2+ 
dose-dependently lowered the activation tempera- 
ture. ***, P < 0.001; n= 8-9. 
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cause an increasing fraction of TRFV1 channels to be 
heat-activated at room temperature. 

We characterized the energetic effect of Mg 2+ on heat 
activation with thermodynamic analysis. It is well estab- 
lished that large changes in enthalpy and entropy are 
hallmarks of the highly temperature-sensitive heat acti- 
vation process of TRPV1 and other thermo-TRP chan- 
nels (Voets et al., 2004; Yang et al., 2010; Yao et al., 
2010). Reflecting a large energetic effect by Mg 2+ on 
heat activation, we found that both enthalpic and entro- 
pic changes were gready reduced when heat activation 
was recorded in the presence of Mg 2+ (Fig. 9 A). At the 
10 mM concentration, Mg 2+ shifted the free energy 
change associated with heat activation by 2.7 kcal/mol. 
Accompanying this energetic effect, there was a sub- 
stantial shift in the temperature activation range at 
higher Mg 2+ concentrations. As shown in Fig. 9 B, the 
half-activation temperature was significandy lowered by 
30 mM Mg 2+ . The large energetic effects on heat activa- 
tion suggest that Mg 2+ may direcdy interact with the 
heat activation process. 

How does Mg 2+ activate TRPV1? 

To understand how Mg 2+ activates TRPV1, we analyzed 
our experimental observations using the multi-allosteric 
model shown in Scheme II (Fig. 1 B). The first step of 
this analysis was to assign the value of equilibrium 
constant L for the C<-»0 transition, which reflects the 
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Figure 8. Shift in TRPV1 heat activation measured by fluores- 
cence imaging. (A) Change in fluorescence intensity as a func- 
tion of temperature in the presence of 1 mM (black) or 10 mM 
(red) Mg 2+ . Arrows indicate the activation threshold temperature. 
(B) Comparison of heat-induced fluorescence increases at two 
Mg 2+ concentrations. Error bars indicate ±SEM. (C) Mg 2+ signifi- 
cantly lowers the activation threshold temperature. ***, P < 0.001; 
n = 58 (for 1 mM Mg 2+ ) and 28 (for 10 mM Mg 2+ ). 



equilibrium state of the channel pore in the absence of 
any activation stimulus. Because spontaneous openings 
of TRPV1 are very rare at room temperature, the equi- 
librium strongly favors the closed state, i.e., L must be 
Indeed, a small L value was taken in all previous 
studies of TRPV1 using allosteric modeling (e.g., 2.6 x 
10 3 in Jara-Oseguera and Mas, 2013; 4.2 x 10~ 4 in 
Matta and Ahern, 2007; 8 x 10" 6 in Yao et al., 2010). In 
the present study, Lwas taken to be 0.005 (corresponding 
to a spontaneous open probability of 0.5%) as an upper 
limit value. A smaller value would better describe the 
data but does not alter the arguments developed here. 

The voltage activation branch 

The equilibrium constant of the voltage-dependent tran- 
sition R<-»A in Scheme II (Fig. 1 B) can be expressed as 



D = exp 



qp(v-v half ) 



RT 



Mg 2+ may potentially affect the R<->A equilibrium by 
changing either the gating charge q or the half-activation 
voltage V ha if. As shown in Fig. 5, q remained unchanged 
in the presence of 100 mM Mg 2+ (0.53 ± 0.03 e 0 vs. 0.56 ± 
0.06 e 0 ). Vhaifi however, shifted substantially to a less de- 
polarized level, from 130.5 ± 1.6 mV to 84.3 ± 1.9 mV. 
The question is whether this near 50-mV shift in V ha ij is 
sufficient to substantially affect the C<->0 equilibrium 
to open the channel. 

The voltage-dependent activation (Fig. 10 A) is de- 
fined as 



LO + JyD) 
l + D + L(l + J V D)' 



(3) 



Using L = 0.005 and the values for D defined by ex- 
perimentally determined q and Vhaifi the value for the 
allosteric coupling factor J v could be determined by fitting 
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Figure 9. Thermodynamic analysis of Mg 2+ effects on heat activa- 
tion. (A) Measured AH values (open bars) and AS values (shaded 
bars) from heat-induced currents. Error bars indicate ±SEM. 
(B) Box-and-whisker plot of the half-activation temperature T hal f 
at the corresponding Mg 2+ concentration. ***, P < 0.001; n = 5-6. 
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the G-V relationship in the absence of Mg 2+ . A rather 
small value of/y= 7 best reproduced the observed open 
probability level and the observed G-V relationship. 
This is anticipated because, in the absence of another 
stimulus, voltage is a poor activator for TRPV1. Even at 
+300 mV, when D was $>l and the R<-»A equilibrium 
fully shifted to the right, the channel open probability, 
P 0 = LJv/ ( 1 + LJv) , remained very low in the absence of 
Mg 2+ (~1.7% at +80 mV; Fig. 5 A, left). Because of the 
small Jv value, a shift of the R<->A equilibrium has a 
minor influence on channel activation, as can be seen 
in Fig. 10 A. In the presence of 100 mM Mg 2+ , the maxi- 
mum open probability predicted by the model is 3.3% 
(at 300 mV), much lower than the experimentally de- 
termined level, which is 27% at +80 mV and room tem- 
perature. Hence, Mg 2+ can't open the channel through 
the voltage-dependent branch of Scheme II (Fig. 1 B) . 
The shift in V hal j observed in the presence of Mg 2+ can be 
explained instead by an allosteric coupling of the R<-»A 
equilibrium to the C<-»0 equilibrium; this will be dis- 
cussed later in "Energetic effects of heat activation on 
the voltage and capsaicin branches." 

The capsaicin activation branch 

As the Mg 2+ effects on channel activation were observed 
in the absence of capsaicin, according to the allosteric 
mechanism of Scheme II (Fig. 1 B) , they can't be achieved 



through the capsaicin activation branch. The observed 
shift in capsaicin dose-response relationship in the 
presence of Mg 2+ is again due to allosteric coupling; this 
will be discussed later in "Energetic effects of heat acti- 
vation on the voltage and capsaicin branches." 

The heat activation branch 

The equilibrium constant C for the heat-dependent 
transition N<-»P is determined by the associated enthal- 
pic and entropic changes: 

C = exp(-A% T + A%). 

A consensus from multiple studies of the thermody- 
namics of TRPV1 gating is that AH and AS are substan- 
tially larger than those observed from most channel 
gating processes (Voets et al., 2004; Yang et al., 2010; 
Yao et al, 2010; Cao et al., 2013). In agreement with 
these previous studies, AH and AS were estimated in 
the present study to be 69.7 ± 4.3 kcal/mol and 217.3 ± 
13.6 cal/mol/K, respectively (Fig. 9 A). The large AH 
and AS values make the N<-»P equilibrium extremely 
temperature-sensitive. When coupled to the C<->0 equi- 
librium, they make the channel activation highly de- 
pendent on temperature. The temperature-dependent 
activation (Fig. 10 B) is defined as 
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Figure 10. Identification of the Mg 2+ target in 
the multi-allosteric TRPV1 gating model. All 
simulations were generated using values derived 
from experimental data presented in Figs. 3, 5, 6, 
and 9, and the following parameters: L = 0.005, 
J H = 300, and J v = 7. (A) Effects of Mg 2+ on channel 
open probability through the voltage-dependent 
branch, (left) The voltage-dependent branch of 
the gating pyramid, with Mg 2+ affecting the equi- 
librium constant D. (right, top) Normalized P-V 
curves in the absence (black curve) or presence 
(red curve) of 100 mM Mg 2+ . Superimposed are 
experimentally determined G-V relationships 
(closed squares), (right, bottom) Predicted chan- 
nel open probabilities using Eq. 3. The maximal 
open probability level in the presence of 100 mM 
Mg 2+ is significantly lower than the experimen- 
tal observation (indicated by the broken line), 
which is ^27% at +80 mV and room temperature 
(see Fig. 3 B). (B) Effects of Mg 2+ on channel 
open probability through the heat-dependent 
branch, (right, top) Predicted heat dependence 
in the absence (black curve) or presence of 10 mM 
(red curve) or 30 mM (green curve) Mg 2+ . Su- 
perimposed are experimentally determined I-T 
relationships (circles), (right, bottom) Predicted 
channel open probabilities using Eq. 4. 
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Because heat can bring the open probability to a high 
level,/// must be 9>l. This can be seen in the reduced 
Eq. 4 when C > 1 (a condition reached at elevated 
temperatures) : 

p= JJh_ 

° 1 + LJh' 

Given that L = 0.005, to reach an open probability of 
50 %,Jh needs to be 200; to reach an open probability of 
90%, J H needs to be 1,800. The large ///value strongly 
couples the C<->0 equilibrium to the N<-»P equilibrium, 
which is expected for a highly temperature-sensitive 
channel. Because of this strong coupling, it is anti- 
cipated that Mg 2+ can effectively activate the channel by 
shifting the N<-»P equilibrium. 

We observed that, in the presence of increasing 
concentrations of Mg 2+ , both the values of AH and AS 
were substantially reduced (Fig. 9 A), which indicates 
a reduction of the overall temperature sensitivity. Also 
noticeable is the substantial shift of the half-activation 
temperature to lower temperatures (Fig. 9 B) , indicat- 
ing that at room temperature the N<-»P equilibrium 
progressively shifts to the right by Mg 2+ . We found that 
both high open probability and left-shifted threshold 
temperature can be satisfactorily reproduced assuming 
Mg 2+ affects only the N<->P equilibrium (Fig. 10 B). This 
apparently is how Mg 2+ activates the channel: by shifting 
the equilibrium of the heat-dependent transition, Mg 2+ 
moves the activation range toward room temperature, 
causing the channel to be heat-activated. 

Energetic effects of heat activation on the voltage 
and capsaicin branches 

The preceding discussions demonstrate that potentia- 
tion of TRPV1 activation by Mg 2+ can be explained by 
a direct effect on the N<-»P equilibrium. To conclude 
the multi-allosteric model analysis, let's consider the 
energetic effects of Mg 2+ on the voltage and capsaicin 
branches of Scheme II (Fig. 1 B) . While Mg 2+ shifts the 
N<->P equilibrium to activate the channel, the resulted 
shift of the C<->0 equilibrium would also make it easier 
for capsaicin and voltage to open the channel. The volt- 
age- and capsaicin-dependent activations in the pres- 
ence of this influence are defined as 

p L{l + J H C + J V D + J H CJ V DJ VH ) ( 5 ) 

° 1 + C + D + J m CD + L(l + J„C + J V D + J H CJ v DJ m ) 



and 



respectively. Fig. 11 shows changes in the G-V curve (A) 
and the capsaicin dose-response curve (B), assuming 
that Mg 2+ only affects the N<-»P equilibrium. Using the 
same set of parameters that generated the results shown 
in Fig. 10, Eqs. 5 and 6 predict that Mg 2+ indeed sub- 
stantially elevates the open probability upon depolariza- 
tion (Fig. 1 1 A) and left-shifts the capsaicin dose-response 
curve (Fig. 11 B). These simulations demonstrate that 
voltage- and capsaicin-dependent activations can be ad- 
equately altered through allosteric coupling of the Mg 2+ 
effects on heat activation. 

Divalent cations induce TRPV1 desensitization to heat 
but not to capsaicin 

Thermodynamic analysis suggests that Mg 2+ activates 
TRPV1 by potentiating heat activation. If this is true, 
Mg 2+ should be able to affect the heat activation process 
independent of the other stimuli. To directiy test this 
prediction, we replaced Mg 2+ with Ba 2+ in our experi- 
ments and examined the channel desensitization pro- 
cess. Like Mg 2+ , Ba 2+ was found to also strongly potentiate 
the channel, eliciting an appreciable current at room 
temperature. As expected, all tested properties of Ba 2+ 
potentiation mirrored that of Mg 2+ (unpublished data) . 
The overall similarity between the effects of Mg 2+ and 
Ba 2+ on TRPV1 indicates that Ba 2+ potentiates the chan- 
nel through the same pathway that mediates the Mg 2+ 
potentiation effect. We noticed, however, that the Ba 2+ - 
elicited macroscopic current was in general much smaller 
in amplitude compared with that elicited by Mg 2+ . The 
observation led us to investigate the desensitization be- 
havior of the channel. 

TRPV1 activity is known to undergo rapid desensitiza- 
tion when recorded in intact cells. The process is Ca 2+ - 
dependent and is mediated by calmodulin (Numazaki 
et al., 2003; Rosenbaum et al., 2004; Lishko et al., 2007; 
Grycova et al., 2008; Lau et al., 2012), phosphatase 
(Docherty et al., 1996; Mohapatra and Nau, 2005), and 
other factors (Lukacs et al., 2007; Ufret-Vincenty et al., 
2011; Cao et al., 2013). Removal of Ca 2+ prevented this 
rapid desensitization process. However, removing Ca 2+ 
revealed a slower, Ca 2+ -independent desensitization pro- 
cess in heat-activated TRPV1 channels, as can be seen 
in the current traces shown in Fig. 7 A. Interestingly, 
Mg 2+ - and Ba 2+ -induced TRPV1 activation also exhibited 
Ca 2+ -independent desensitization. Hints of a slow de- 
sensitization process can be seen in the representative 
fluorescence recording shown in Fig. 2 B and current 
recording shown in Fig. 3 A. Ba 2+ -induced channel de- 
sensitization was much more rapid than Mg 2+ -induced 
desensitization. An example demonstrating the time 
course of Ba 2+ -eli cited current is shown in Fig. 12 A, 
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where a rapid decline in current amplitude is seen with 
a dominant time constant of 5.8 ± 1.0 s (n = 4). The 
Mg 2+ -elicited current declined with a dominant time 
constant of 84.5 ± 17.9 s (n = 3), more than 10 times 
larger than that of Ba 2+ -induced desensitization. It is 
hence most likely that rapid desensitization substan- 
tially diminished the Ba 2+ -elicited current, making the 
Ba 2+ potentiating effect appeared to be less strong than 
that of Mg 2+ . 

To investigate the underlying mechanism for the di- 
valent cation-induced Ca 2+ -independent desensitization 
process, we designed a two-step protocol, as shown in 
Fig. 12 B. A short temperature ramp was first delivered 
to induce heat activation in cell-free inside-out patches 
in a Ca 2+ -free solution. After closing the channels by 
cooling, the same membrane patch was challenged by 
3 uM capsaicin to induce agonist activation. In the ab- 
sence of divalent cation, temperature-induced current 
kept increasing in amplitude up to 50°C, yielding a 
current comparable in size to that induced by capsaicin 
(Fig. 12, B and D). However, for patches exposed to 
Ba 2+ before and during recording, heat could no lon- 
ger effectively activate the channels, though the capsa- 
icin-induced current remained (Fig. 12, C and D). The 
lack of heat-induced current suggests that during the 
extended presence of Ba 2+ , the heat activation machin- 
ery was mostiy desensitized. Indeed, as can be seen in 
Fig. 7 A, desensitization of heat-activated channel cur- 
rent becomes more prominent and occurs at progres- 
sively lower temperatures in the presence of increasing 
concentrations of Mg 2+ . Desensitization of the heat 



activation machinery, however, did not prevent capsaicin- 
induced activation (Fig. 12 E), which is consistent with 
the notion that capsaicin activation and heat activation 
use separate pathways. These results further confirmed 
that Mg 2+ and Ba 2+ potentiate TRPV1 by selectively in- 
teracting with the heat activation pathway. 

DISCUSSION 

Our data demonstrate that the divalent cations Mg 2+ 
and Ba 2+ dose-dependently potentiate TRPV1 activa- 
tion. This gating effect is partially covered by an inhibi- 
tory effect on the single-channel conductance. Adjusting 
for the conductance difference revealed strong poten- 
tiation of TRPV1 gating by Mg 2+ that brings the channel 
open probability to ^25% at room temperature. The 
substantial effect on channel gating, observed from 
both current and fluorescence recordings, are in close 
agreement with several previous studies using divalent 
and multivalent cations (Ahem et al., 2005, 2006; 
Tousova et al., 2005; Riera et al., 2007; Ohta et al., 2008; 
Luebbert et al., 2010), with minor quantitative differ- 
ences that are accountable to variations in experimen- 
tal conditions. For example, the apparently lower EC50 
value for Mg 2+ in a previous study (Riera et al., 2007) 
might be due to a combination of differences in chan- 
nel homologues (mouse vs. human), Ca 2+ indicator dye 
(Fluo-4 vs. Fura-2), and perhaps the use of SC>4 2 ~ in- 
stead of CP as the counter ion in the previous study; 
indeed, at the same concentration (10 mM), NiSC^was 
found to be more effective than NiCl2 in potentiating 




1.0 



0.8 



0.6- 



0.4 



0.2 




30 mM 



0.0-1 1 1 =1 r 

•200 -100 0 100 200 300 
Voltage (mV) 




100 1000 
[Capsaicin] (nM) 



Figure 1 1 . Mg 2+ potentiation of heat activation can 
substantially affect voltage- and capsaicin-dependent 
activation. All simulations were generated using val- 
ues derived from experimental data presented in 
Figs. 3, 5, 6, and 9, and the following parameters: 
L = 0.005, J H = 300, J v =7,Jc= 800, J m = 1, and J m =l. 
(A) Mg 2+ effects on voltage-dependent activation, 
(left) The dual-allosteric model containing the 
heat and voltage branches, with Mg 2+ affecting the 
equilibrium constant C. (right) Predicted voltage- 
dependent activation using Eq. 5. The broken line 
indicates the position of +80 mV. (B) Mg 2+ effects 
on capsaicin-dependent activation, (left) The dual- 
allosteric model containing the heat and capsaicin 
branches, with Mg 2+ affecting the equilibrium con- 
stant C. (right) Predicted capsaicin dose-response 
curves using Eq. 6. 
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TRPV1 (Luebbert et al., 2010). More importantly, our 
present study suggests that potentiation of TRPV1 by 
Mg 2+ and Ba 2+ is achieved predominantly by lowing the 
channel's heat activation threshold so that the channel 
starts to be heat-activated at room temperature. 

Although phenomenally divalent cations shift the 
temperature range at which heat activation occurs, 
quantitative comparison of the energetic effects of Mg 2+ 
on each activation pathway suggests that mechanisti- 
cally divalent cations potentiate TRPV1 by reducing the 
energy barrier for heat activation. Both the enthalpic 
change and the entropic change associated with heat 
activation are substantially lowered in the presence of 
Mg 2+ . As a result, it requires less thermal energy to drive 
heat activation. Thermodynamic studies showed that 
heat activation of thermo-TRP channels is associated 
with large changes in enthalpy and entropy (Voets et al., 
2004; Yang etal., 2010; Yao etal., 2010). The magnitude 
of these changes indicates a substantial conformational 
change in the channel protein (Clapham and Miller, 
201 1) . Although the nature of this heat-induced confor- 
mational change remains elusive, its presence could be 
detected directly with fluorophores attached to the 
pore turret (Yang et al., 2010). Interestingly, while fluo- 
rophores on the pore turret reported a conformational 
change during heat-dependent activation, there is no 
evidence of a conformational change at the same posi- 
tion during capsaicin- or voltage-dependent activation. 
The large change in thermodynamic properties by Mg 2+ 
observed in the present study is consistent with these ear- 
lier findings, suggesting that the presence of Mg 2+ and 
Ba 2+ may promote the heat activation conformational 



change through which it exerts potentiation effects on 
TRPV1 activation. It is important to note that changes 
in thermodynamic properties may arise from a change 
in either the equilibrium constant C of the N<->P transi- 
tion (which may be considered to represent the "heat 
sensor") or the coupling factor J u . Measurements in the 
present study could not distinguish between these two 
alternatives. Future experiments are required to estab- 
lish whether Mg 2+ directly affects movement of the heat 
sensor or the structure that couples the heat sensor to 
the activation gate. For this reason, in this paper we 
restrict our discussion to distinguishing between heat, 
voltage, and capsaicin "activation pathways." 

TRPV1 is a polymodal cellular sensor for voltage, cap- 
saicin, and heat (Tominaga et al., 1998) . These stimuli 
likely work through separate activation pathways (Jordt 
et al., 2000; Matta and Ahern, 2007; Grandl et al., 2010; 
Yang et al., 2010; Cui et al., 2012) but synergistically ac- 
tivate the channel. As a result of functional coupling, 
experimental perturbation to one modality would al- 
losterically affect the channel's response to other mo- 
dalities. While mutagenesis studies in recent years have 
yielded several important observations that may shed 
new light on the structural basis for heat activation, it 
is important to realize that, because of functional cou- 
plings between different activation pathways, almost 
every perturbation to channel gating may manifest in 
a change of the channel's response to other stimuli 
(Brauchi et al., 2004; Latorre et al., 2007). Current po- 
tentiation by divalent cations demonstrates again the 
synergistic nature of activating a polymodal channel 
like TRPV1 . It underlines the need to distinguish direct 
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Figure 12. Divalent cations induce 
desensitization of the heat activation 
pathway but not the capsaicin activa- 
tion pathway. (A) Normalized current 
traces demonstrating transient activa- 
tion of TRPV1 by 130 mM extracellular 
Ba 2+ (left) or Mg 2+ (right), followed by 
rapid desensitization. (B) Time course 
of whole-cell current (black trace) at 
+80 mV evoked first by heat (red trace 
below current trace) and subsequently 
by 3 pM capsaicin (black bar above cur- 
rent trace). (C) Same as for B, but in 
the presence of 130 mM Ba 2+ instead 
ofNa + . (D) Box-and-whisker plot of the 
ratio between current evoked by rais- 
ing temperature to 45°C and current 
evoked by 3 pM capsaicin in a bath so- 
lution containing the indicated cations. 
***, P < 0.001; n = 4-5. (E) Schematic 
gating diagram illustrating the separa- 
tion of the capsaicin activation pathway 
and the heat-Mg 2+ -Ba 2+ pathway (left). 
In the presence of Mg 2+ or Ba 2+ , the 
heat-Mg 2+ -Ba 2+ pathway is desensitized 
over time but the capsaicin activation 
remains functional (right). 
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versus indirect interactions and highlights the advantage 
of modality-specific perturbation methods in mecha- 
nistic investigations. 

Desensitization of the TRPV1 gating machinery 
TRPV1 exhibits complex desensitization properties. The 
dominant and rapid desensitization process is found 
to be Ca 2+ -dependent (Docherty et al., 1996; Koplas 
et al., 1997). It is shown that Ca 2+ -dependent desensiti- 
zation is partly caused by channel modulation by Ca 2+ - 
calmodulin, which binds to two intracellular sites, one 
in the channel's C-terminal region and the other in the 
Ankyrin-like repeat domains of the N-terminal region 
(Numazaki et al., 2003; Rosenbaum et al., 2004; Lishko 
et al., 2007; Grycova et al., 2008; Lau et al., 2012). 
Desensitization occurs when Ca 2+ enters the cell through 
activated TRPV1 channels and activates calmodulin. 
This Ca 2+ -dependent desensitization process is often 
observed in experiments using whole-cell recordings 
but is also observed in cell-free patch recordings when 
Ca 2+ is present. It plays a major role in channel down- 
regulation upon prolonged or repetitive stimuli under 
physiological or pathological conditions (for review see 
Vyklicky et al., 2008). 

Preventing Ca 2+ -dependent desensitization by exclud- 
ing Ca 2+ from experimental solutions reveals the exis- 
tence of other types of desensitization. For example, 
repetitive or extended heating to high temperatures 
leads to current decline in TRFVT (Tominaga et al., 
1998) , which is demonstrated in Fig. 7 A of this paper. We 
observed a similar Ca 2+ -independent desensitization pro- 
cess in Mg 2+ -and Ba 2+ -induced currents. Because Mg 2+ 
and Ba 2+ shift the heat activation process, it is likely that 
desensitization of TRPV1 after heat or divalent cation ac- 
tivation reflects a common gating process. This is con- 
firmed by the observation that Mg 2+ /Ba 2+ -desensitized 
channels can no longer be heat-activated. The finding 
that capsaicin remains effective in activating TRFV1 after 
Mg 2+ - or Ba 2+ -induced desensitization is fully consistent 
with this idea and confirms the separation of the capsa- 
icin and heat activation pathways. More importantiy, the 
observation suggests that heat- and divalent cation- 
induced desensitization does not involve the closure of 
a distinct inactivation gate. Instead, it is an earlier step of 
the heat activation process that has become desensitized 
to heat and divalent cations so that it can no longer sus- 
tain the opening of the activation gate. 

Several molecular events are known in various channel 
types to terminate ion flux upon extended activation. 
Two inactivation processes occur in Shaker potassium 
channels: the N-type inactivation involves physically plug- 
ging the open pore with an "inactivation ball" peptide 
(Hoshi et al., 1990; Zagotta et al., 1990), whereas the 
C-type inactivation arises from conformational changes 
in the outer mouth of the pore (Choi et al., 1991; Yellen 
etal., 1994; Baukrowitz and Yellen, 1995, 1996; Liu etal., 



1996). It is noted that for both inactivation processes a 
distinct "gate" is shut to prevent ion flux even though the 
main activation gate remains open. Ligand-gated chan- 
nels terminate ion flux in the extended presence of li- 
gand molecules through a desensitization process that 
recloses the same gate that is opened during activation. 
In the glutamate receptor this is achieved by relaxing a 
part of the gating machinery that experiences a strain 
caused by the activation stimulus (Sun et al., 2002). Hy- 
perpolarization-activated cyclic nucleotide-modulated 
(HCN) channels desensitize to voltage through a slip- 
page in the coupling between the voltage sensor and the 
activation gate that relieves the strain in the gating ma- 
chinery (Shin et al., 2004). An important distinction be- 
tween the inactivation mechanism and the desensitization 
mechanism is that an inactivated channel cannot reopen 
before the inactivation gate opens again, whereas a de- 
sensitized channel can be opened by a stimulus that 
avoids the desensitized gating machinery and works 
through a separate pathway. In the present study we pres- 
ent direct evidence in support of this prediction. Because 
Ca 2+ -independent desensitization of TRPV1 does not 
prevent subsequent capsaicin activation, it is most likely 
not caused by a distinct inactivation gate, but rather by 
closing the main activation gate that is formed by the 
lower part of the S6 segments (Salazar et al., 2009). The 
desensitization process may occur locally to the same 
channel structure that undergoes heat-induced activa- 
tion conformational changes, leading to the closure of 
the activation gate. 
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